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Abstract 

We examine the effect of the intra- and interspecies scattering lengths on the dynamics of a 
two-component Bose-Einstein condensate, particularly focusing on the existence and stability of 
solitonic excitations. For each type of possible soliton pairs stability ranges are presented in 
tabulated form. We also compare the numerically established stability of bright-bright, bright- 
dark and dark-dark solitons with our analytical prediction and with that of Painleve- analysis of 
the dynamical equation. We demonstrate that tuning the inter-species scattering length away from 
the predicted value (keeping the intra-species coupling fixed) breaks the stability of the soliton pairs. 

PACS numbers: 03.75.Kk, 03.75.Mn 



* apagyi@phy.bmc.hu 



1 



I. INTRODUCTION 



Since Bose-Einstein condensates (BECs) can be routinely prepared in laboratories, ul- 
tracold gases became a very important testbed for many predictions of condensed matter 
physics The experimental examination of binary condensates started nearly the same 
time as for single condensates by using two different quantum states of the same species, such 
as 87 Rb jsj or 23 Na 3|. With the development of sympathetic cooling ultracold mixtures have 
been assembled from two different alkalies, 41 K- 87 Rb {4,5)], 7 Li- 133 Cs and 87 Rb- 133 Cs [l| 
or for different isotopes of the same alkali atom 85 Rb- 87 Rb The tunability of the inter- 
and intra-species scattering lengths via drivingthe mixture through a Feshbach-resonance 
has also been experimentally demonstrated [JO]. 

The ability to create Bose-Einstein condensate (s), a highly coherent form of matter, also 
facilitated the convergence of two fields of physics: condensed matter physics and quantum 
optics, and therefore BECs became favourable candidates for examining the effects of nonlin- 
earity in matter waves, where this nonlinearity originates from the mean-field representation 
of the interatomic interaction. The similarity between electromagnetic waves in nonlinear 
medium and coherent matter waves is also expressed in the equations of motion which is the 
nonlinear Schrodinger equation (NLS) for the former and the Gross-Pitaevskii (GP) equa- 
tion for the latter. Although the physical interpretation of these equations is different, their 
structures are the same apart from the external potential term. Furthermore, in some cases 
this extra term can even be removed \^ and the GP equation is transformed into a form 
coinciding with the NLS equation. Consequently all results for the NLS equation known in 
nonlinear optics can be readily adapted to Bose-Einstein condensates. 

One of the surprising phenomena of nonlinear optics is the existence of particle-like wave- 
forms, the so-called solitons 10| . Such excitations have already been experimentally observed 
e- or two-component Bose-Einstein condensates: dark solitons flll.[l2|. bright solitons 



m sing 



13 



14], their two-component coupled analogues, the dark-dark [15J, bright-bright |16| or 



even dark-bright 17], ll8[ multi-component solitary waves 18 



19. 



However, the question of existence of solitons needs more attention than simply recog- 
nising the similarity between the two governing equations. The existence of solitons is 
strongly related to the integrability of the given physical model. An usual test to deter- 
mine whether or not an equation is integrable is the Painleve-test (P-test) 20j. It was 



shown that without any potential term or inhomogeneity the one- dimensional NLS equa- 



tion, iu t + u xx ± 2 \u\ u 
or multi-component 22 



0, is completely integrable both in its one-component 
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2l| 



23[ form. However, the inclusion of a potential term, v(x,t)u, in 



the one-component NLS equation or different coupling strengths in the multi-component 
NLS/GP equations fundamentally changes their integrability 2% l25|. It was shown that 



integrability is preserved provided the external potential, v(x,t), has a specific form 25 |. 



The authors examined 26[ the integrability of the two-component coupled Gross-Pitaevskii 
(CGP) equations and lead to similar conclusion: the scattering lengths and the external 
potentials cannot be arbitrary if the integrability of CGP is to be preserved. The inter- and 
intraspecies scattering lengths must satisfy the following equation 

266-^i6~^6 _ Q + 1) 2 + 7 

66 - «i« 2 16 

where 6 = 011/021, 6 = 022/012, «i = aWa«21, ^2 = ^22/^12 and /i^ denotes the reduced 
mass of pair of particles composed by an atom from the ith and j'th species. On the right 
hand side of Eq. (pQ) n is a non-negative integer number. One may call n a classification 
number, because it determines the form of the external potentials for which CGP equations 
remain integrable. For example, for n = 2, the external potential, apart from the quadratic 



trapping potential, may even contain an imaginary time- dependent term [26(. This term 
can mimic the loss or gain in the number of particles of the given species. We note here that 
usually dissipation works against long-living coherent matter waves, however, the importance 



of this imaginary potential term has been analysed in 
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soliton solution 



and shown to permit exact 



]. In the context of BECs at finite temperature, the interaction of the 
condensate with the thermal cloud could also be taken into account as an imaginary term 
in the governing GP equation. This interaction, due to its stochastic nature, can influence 
the dynamics of the solitons, via density-fluctuation. 

In this paper we carry out an analysis on how the intra- and interspecies interaction 
influences the dynamics of a binary mixture of Bose-Einstein condensates. We select out of 
the many possible systems of two component BECs the pairs 87 Rb- 87 Rb (prepared in two 
distinct hyperfine states), 23 Na- 87 Rb, and 7 Li- 39 K. In the former two systems there is a 
possibility to study stability of bright-dark and dark-dark soliton pairs, while the last pair 
is capable to sustain bright-bright and bright-dark excitations. 

The organization of the paper is as follows. Section [Til defines a quasi one-dimensional 



model derived from the general three-dimensional, coupled Gross-Pitaevskii equations as- 
suming cigar-like harmonic oscillator trap potential. In the first part of section IIHI we shall 
perform a stability analysis based on coupled soliton excitations, and on Eq. ([I]) of the 
P-test. In the second part of section IIII[ possible new interesting modes exhibited by the 
bright-bright solitons will be shown. Section [IV] is devoted to a conclusion and the summary. 



II. THE MODEL 



In ultracold gases the interaction between two particles can usually be well described by a 
scalar parameter, the scattering length. For a two-component Bose-Einstein condensate one 
has to introduce three, possibly different, scattering lengths, characterising the intra-species 
interactions (an, 022) and the inter-species coupling (au = 021)- 

In the mean-field approximation a two-component BEC is described by the coupled Gross- 



Pitaevskii equations 
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3l| in the form: 

h 2 



d 



2 m,- 



A 



Vi 



(2) 



where m, denotes the individual mass of the ith atomic species, fiy = 2nh 2 aij/ fiij with 
being the 3D scattering length, = mimj/(mi + rrij) is the reduced mass, and Vi denotes 
the external trapping potential. In the following, indices % and j label the components, 
therefore, take only two values, 1 and 2. In the case of real trap potentials the normalisation 
of the wave functions reads as iVj = Jl^ifdV with Ni denoting the number of atoms in the 
ith component. We exclude those cases from our analysis where the species can transform 
into each other, therefore the number of atoms in each component hereafter is conserved. 



A. Transformed equations 



If the three dimensional quadratic trapping potential is weak in one direction, i.e. 



Vi = -rrii 
1 2 



2 2, 2 / 2 1 2\ 

"i,x x + u i,± (y + z ) 



(3) 



where Ui X Ui±, one may replace the three-dimensional equations (j2J) with a coupled system 
of quasi one-dimensional GP equations. Although Eq. (El) is a nonlinear equation, physically 



we may assume that the weak x-direction decouples from the strong [yz\— plane, therefore 
the macroscopic wave functions can be written as 



(4) 



with Xt,± represents the ground-state solution of the corresponding two-dimensional Schrodinger 
equation in the [yzj-plane. The external potential introduces suitable units of length and 



time as a± = \7h/m\ui\^_ and r = l/u^x, respectively. By rescaling the spatial and 
temporal variable with a±_ and r one obtains two quasi one-dimensional GP equations 



#2,i 



M 

2 



^d xx + ^x 2 + b ll \^j 1 \ 2 + b 12 \^j 2 \ 2 



1 
2 

2 2k 



621 IV'll 2 + ^22 |"02| 2 



^1, 
^2, 



(5a) 
(5b) 



where b u = 2anNi, b 22 = 2a 22 Nin/'j, fo 12 = 621 = 2a 12 iVi(l + «)/(l + 7), 7 = W2,±A>i,±, 
k = m 1 /m 2 , Ai = Ux tX /uJi t j_, A 2 = W2,as/ W i,±- The normalisation is such that J |^i| 2 <ir = 1 
and J" 1^2 1 2 rfa; = N 2 /Ni. Moreover, the relation = k must hold if both species experience 
the same harmonic potential. (Note the slight departure from Ref. 32[ in the definition 
of 622 which, however, may result in large difference of values of 622 if a two-component 
condensate contains species with different masses mi 7^ m 2 .) 



B. Thomas-Fermi background 

If the kinetic energy term is negligible compared to the potential energy terms in (I5al b) . 
then one may apply the Thomas-Fermi approximation to determine the density distribution 



of the ground state. Following |33[ we write the corresponding wave-functions as 



resulting in TF densities 



and TF energies 



if>i(x, t) « $J F (x) exp (-iEj F t/h), (6) 



E? F = (bnAtxl + b 12 A 2 x 2 2 ) /A, (8a) 
^ 2 TF = {buA^ 2 + b 22 A 2 xj) /A, (8b) 
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where A = 611622 — b\ 2 . The parameters Ai and X{ represent the amplitude of the density 
and the extension of the condensates, respectively. All these quantities are determined by 
the system parameters bij, N iy Aj according to the following relations: 

A 6 22 . 2 h2 > 2 , , bn 2 b 12 2 / n v 

1 = T 1 ~ 2k" 2 ' = 2k" 2 ~ T x ' ™ 

while the extensions are 

xi= (it) and x2= {ItM) • (io) 

Although the Thomas- Fermi density distribution is not physical at x — Xi, it still provides 
a good starting point for analytical calculations. In our numerical treatment we will not use 
this approximation, rather start our simulations from the appropriate ground state solution 
of the one-dimensional GP equation. 



C. Coupled soliton excitations 



Now we are seeking solutions of Eqs . fl5al -b) which support soliton excitations. A static 
soliton excitation can be written as [33] 



*l>i(x,t) = $7 F (0)v9,(x)exp(^t). 



(11) 



By inserting this ansatz into the GP Eqs. (JSaj-b) and neglecting the small potential contri- 
butions one obtains the coupled soliton equations as follows 



-- d xx + 6n 1 9^1 1 2 + 612 IH 2 



E 2 (p 2 = ~^d xx + b 21 \<fii\ 2 + b 22 |y? 2 | 2 



<P2 



(12a) 
(12b) 



with b^ = bijAjX^/A. The normalisation of the soliton solutions reads as follows 



\(px\ dx 



\(p 2 \ dx 



A 



A,xV 

A N 2 
A 2 x\ Ni ' 



(13a) 
(13b) 



-L 2 ^2-^2 

where the integrations, in both cases, are over the spatial extension, L iy of the solitons. 

The above coupled equations admit generic moving soliton solutions of the types: bright- 
bright (BB), bright-dark (BD), dark-bright (DB) and dark-dark (DD). We shall investigate 



(3 



here a simple static bright-dark soliton pair solution by taking the first component to be a 
static bright soliton 



ip BD (x) = q 1 sech (hx), tp BD (x ->■ ±00) = 0, (14a) 

and the second component to be the static dark soliton 

<pf D (x) = q 2 tanh (k 2 x), <f 2 D (x -»■ ±00) = ±q 2 (14b) 

with yet unknown complex amplitudes g« and wave- vector, fcj. The latter one is related to 
the width of the soliton fcj ~ I /Li. By inserting the above ansatz into the Eqs. fll2al b) 
and equating the coefficients of the constant and ^-dependent terms, respectively, one may 
conclude that the wave-vectors of the dark and bright solitons must be equal, k\ — k 2 = k. 
The amplitudes are expressed by the system parameters as 

k 2 

\<li\ 2 = -; — 0(^12 - 622), (15a) 
Aix\ 



1 2 _ ^ 2 
" A 2 x\ 



l&l = ~rz ^(Kb n - b 12 ). (15b) 
The energy of these excitations read as 



gee = k *(^hi bl2 _ , (16a) 

771BD ; 2 — ^12 , /, Gl x 

#2 = k £ — 6 2 2- (16b) 

Suppose now that the width parameter k is a real number. The modulus of the amplitudes 
of the bright-dark soliton pair superimposed on the Thomas- Fermi background must be real 
numbers, thus one obtains a the following set of conditions for the existence of this bright- 
dark soliton excitations 

Cl = K6l2 A &22 > and C 2 = Khl ~ 612 > 0. (17) 

One may apply the same method to generate static bright-bright or dark-dark soliton 
excitations. In the bright-bright case one obtains the following solutions 

£ 

(fi B (x) = —== — y nb 12 - b 22 sech (kx), (18a) 
k 



(f2 B (x) = - y/b 12 - Kb n sech (hx), (18b) 

\JA 2 x 2 



with the conditions 

d > and C 2 < 0, (19) 
while the energies are Ef B = E^/k = —k 2 /2. The dark-dark coupled soliton solutions read 



as follows 



k 



r A[x 1 



— K ^i2 tanh (kx), (20a) 



k 

(f2 B (x) = _ y/nhn - 6 U tanh (fcx), (20b) 
V^2 ^2 



with the conditions 



d < and C 2 > 0, (21) 

and energies E® D = E^/k = k 2 . Our interesting result shows that in the bright-bright and 
dark-dark cases the energies are uniquely determined by the wave- vector and the mass ratio. 
Note that the existence conditions ( |T71) . (fl9~i) . and (12ip are just the same as obtained in Ref. 
321 ] for the existence of moving soliton pairs, while the constraints for static excitations were 
published in [33] . 



III. STABILITY TESTS BY SIMULATION 



Below we are going to numerically investigate the stability of soliton pairs. To solve 
the time dependent coupled Gross-Pit aevskii equations fl5al -b). a third-order accurate split- 
step Fourier transform method is used as described in Ref. for a single-component 
condensate. Here we solve the time independent coupled Gross-Pitaevskii equations for 



their numerically exact ground states using imaginary time method 35[] combined with the 
split-step operator technique. Choosing initial distributions is necessary to this method, 
and the Thomas-Fermi approximate solution proved to be an effective initial guess for this 
purpose. 

The procedure explained in the previous section can be generalised for solitons moving 
with velocity v. Such a solution is given by 



ft 

V>2 



D = $^ F (0)^ D i 

BD 



X 



vt) exp ( -iEf D t ) exp 



vH 



V(X 



Vt) 



'.VC 2 

I V 



BD/ 



X X 



Vt) 



exp(-^ 2 BD t)exp(-*^^, 2 t 



,(22a) 
(22b) 
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At v = we obtain the static soliton excitation solution as given by ( ITT]) in the preceding 
section. 

The existence conditions (1171) prescribe various relations between domains of the inter- 
and intra coupling strengths, bij, where it is possible to create bright-dark soliton pairs. 
These domains are listed in Table [TJ Although the creation of a bright soliton is generally 
associated with attractive interaction (pa < 0) among the particles, Table [TJ clearly shows 
that, due to the appropriate other couplings, it may also be possible to create a bright soliton 
in case of repulsive interaction (bu > 0). Such a situation occurs in case 4a of Tab. [TJ which 
we are going to analyse below. 



A. Moving bright-dark soliton pairs 

Let us now investigate the stability of a bright-dark soliton pair for the case 4a of Table [TJ 
by considering two experimentally accessible two-component BECs. The first is composed of 
the two hyper-fine states of 87 Rb atoms, the second is obtained from 23 Na and 87 Rb atoms. In 
both cases the scattering lengths are well known and can be tuned over a broad limit by using 
the Feshbach resonance method. Our aim here is to explore the sensitivity of the temporal 
evolution of soliton pair when the inter-atomic coupling strength, 6 12 , is varied around the 
value prescribed by the ratio formula ([1]) obtained by performing a Painleve-analysis of the 



coupled GP equations 26]. In this respect we fix the intra-atomic strengths, ba, to values 
easily accessible to the experiments and vary the inter- atomic interaction values, bu, within 
a small range allowed by the case 4a listed in Table [TJ In order to represent a more realistic 
situation, a small velocity of v — 0.04 is given the solitons, a weak harmonic trapping 
potential is added along the longitudinal direction and the solitons are superimposed on the 
ground state density distribution. This procedure spatially confines the solitons without 



affecting their essential stability properties [32 ] . 



In Fig. [TJ we plotted a bright-dark soliton pair composed of two hyperfine states of 
87 Rb atoms. The top two panels show the case when the inter-species interaction is chosen 
according to Eq. (TTJ. Both solitons oscillate in the harmonic trap in a stable manner with an 
angular frequency slightly less than u x / what the one-dimensional Thomas-Fermi model 



predicts for a single dark soliton 361 ] . The difference is probably caused by the presence 



of the bright soliton, since the bright component fills the dip of the dark soliton, therefore 



Case 



hi b 2 2 



Constraint on 612 



1 + — no bright-dark soliton pair 

2 - + bnK < < 

3a - - -\/6n&22 < &12 < hiK if « 2 &n > ^22 
3b bun < < -\/&n&22 if « 2 &n < ^22 

4a + + V&11&22 < &12 < ^22/^ if k 2 6h < 622 
4b 622/K < 612 < \/&ii&22 if « 2 &n > ^22 

TABLE I. Various domains of inter- and intra-atomic interaction strengths, hj, permitting the 

existence of a bright-dark soliton pair. 
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FIG. 1. (Color online) Evolution of the bright-dark solitons of two-component BEC composed of 
two hyperfine states of the 87 Rb atom. Intraspecies scattering length a\2 = 5.5 nm (612 = 13.6) 
satisfying ratio formula ([1]) (top two panels), a\2 = 5.4 nm (612 = 13.35) disobeying the ratio 
formula (pQ) (bottom two panels). Other parameters an = 0.8 x 5.5 nm (bu = 10.88), 022 = 
1.2 x 5.5 nm (6 2 2 = 16.32), m 1 = m 2 = 87 au, u h ± = 2vr x 710 Hz, Ai = A 2 = 0.2, v = 0.04, 
Ni = 500, N2 = 6600. The snapshots are depicting the soliton pairs at t = 20 and t = 200, 
respectively. 



;he dark soliton has to drag this extra mass as well. This effect has recently been observed 
181 ] with a 87 Rb- 87 Rb condensate, prepared in the \F = 2, nip = 0) (bright soliton) and 

\F = l,mp — 0) (dark soliton) hyperfine states. 

However, if a\2 (or equivalently 612) is tuned away from this particular value, the stability 

is lost, and the initial forms of the solitons are destroyed by the destructive interference of 
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FIG. 2. (Color online) Evolution of a bright-dark soliton pair in a two-component BEC composed of 
23 Na and 87 Rb atoms. Intraspecies scattering length a±2 = 6.831 nm (612 = 7.25) satisfies formula 
(HJ) (top two panels), while a±2 = 3.769 nm (612 = 4) deviates from Eq. (P) (bottom two panels). 
Other parameters are: an = 2.7 nm (pu = 3.4), 022 = 5.5 nm (622 = 3.6), mi = 23 au, 777,2 = 87 au, 
lo 1i± = 2vr x 710 Hz, Ai = A 2 = 0.25, v = 0.15, iV"i = 500, N 2 = 58000. The snapshots are taken at 
t = 20, and t = 2, respectively. 

the constantly emitted and re-captured sound waves. It is worthwhile to mention, although 
only as a qualitative statement, that the appearance of sound waves made the solitons' 
oscillation faster (see the different range of time in the top two and bottom two panels of 
Fig. [p. The sound waves travel faster than the solitons, and after being reflected back 
from the edge of the condensate, they collide with the solitons. The subsequent collisions 
possibly speed up the oscillation and turn it into an irregular sloshing. Despite the irregular 
movement of the solitons, the dark component still captures the bright soliton during the 
motion. 

In Fig. [2] the evolution of a bright-dark soliton pair in binary BEC composed of 23 Na and 
87 Rb atoms is plotted with corresponding snapshots of the density distribution. Upper two 
panels exhibit the bright ( 23 Na) and dark ( 87 Rb) excitations, respectively, which are stable 
for a long period due to the fine tuning of scattering lengths (i, j = 1, 2) which satisfy the 
P-test formula (P). On the lower two panels of Fig. |2]we see, however, that the initial soliton 
excitations do not remain stable but are rapidly dissolved due to detuning the interspecies 
scattering length from the value obeying the integrability condition expressed by Eq. 

dTj- 
ii 



Summarising, both examples exhibited in Figs. [T] and |2] show that a long lived stability 
of the bright-dark soliton pairs can be achieved only in case if the interatomic interaction 
parameters (or 6 12 ) is chosen according to Eq. ([I]). This observation emphasises that there 
may be situations in the creation of two-component BECs when fine tuning of scattering 
lengths according to the P-test formula ([I]) may prove useful. 

At the end of this section we briefly mention an early investigation of the stability of 
a heterogeneous two-component Bose- Einstein condensate by Law et al. {stJ. This work 
carried out a linear stability analysis by calculating the lowest eigenvalue of an excitation. In 
this approach the appearance of a negative eigenvalue signals instability. It was shown for a 
sodium-rubidium system, just as above, that stability occurs only in finite range of the inter- 
species scattering length. The direct quantitative comparison with our work, however, is less 
straightforward because Law et al.'s analysis assumes a spherically symmetric condensate, 
while we analyse a quasi one-dimensional model. 

B. Static bright-bright soliton pairs 

Let us now investigate the temporal stability of a static (v = 0) bright-bright soliton pair 
obtained as the exact solution of Eqs. (l5al -b) in the absence of a trapping potential (Aj = 0). 
The solutions have the form as in Eq. ffTTj) where (fi(x) are chosen to have the functional 
forms given in Eqs. (I18al b). 

The existence conditions (fT9l) prescribe the relations between domains of the inter- and 
intra-species coupling strengths, 6^, which are listed in Table [TTJ Moreover, the common 
wave- vector is k = l/(2Ci). 
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TABLE II. Various domains of inter- and intra-atomic interaction strengths, 6y, permitting the 
existence of a bright-bright soliton pair. 
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FIG. 3. (Color online) Evolution of a bright-bright soliton pair in a binary BEC composed of Li 
and 39 K atoms. Intraspecies scattering lengths and atom numbers are a\2 = 0.2 nm (612 = 0.46), 
N2 = 2029 (first two panels), a\i = 0.275 nm (612 = 0.64) N2 = 2270 (second pair of the panels), 
012 = 0.3 nm (612 = 0.7) N2 = 2347 (third pair of the panels). Other parameters are: a\\ = — 1.4 nm 
(611 = —3.93), (Z22 = — 0.9nm (622 = —1.07), mi = 7au, 1712 = 39au, u)\ t ±_ = 2ir x 710Hz, 
^1 = A2 = 0, v = 0, Ni = 2000. The snapshots are taken at t = 35, t = 35, and t = 20, 
respectively. 
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Scenario 4, for example, describes two condensates for which the intra-species interactions 
are repulsive. This situation, using Hartree-Fock calculation, has been theoretically analysed 
soon after the observation of overlapping condensates prepared from two hyperfme states 



of 87 Rb [2|. For this case, i.e. mi = m 2 , it was established that the two condensates cannot 
co-exist if | a-12 1 < v /a n a 22- Our approach reproduces and extends this result for the case 
of different species. This surprisingly simple relation can be understood using energetic 
arguments; if 612 overcomes the geometric mean of the intra-species interaction strengths, 
the repulsion between the two condensates will separate the two condensate completely and 
they will not overlap any more. However, if the two species attract each other enough, i.e. 



b\2 < — V^n^22, the attraction will dominate and can counteract the individual repulsion 
present in each component. 

Another interesting scenario here is the one listed under the case 3a in Table [Til showing 



that it is possible to create a bright-bright pair, within the range < b\ 2 < in spite 

of the repulsive inter-atomic interaction. In order to investigate the stability of the bright- 
bright soliton pair in this domain we simulate the temporal evolution of the BEC system 
composed of 7 Li and 39 K atoms accessible for experiments. As Fig. [3] shows, depending on 
the value of the interspecies interaction, two types of instability may occur. At values of 
a 12 less than a critical value of about 0.275 nm, the two standing solitons begin to repel 
each other (first two panels) and depart from each other as they were particles, preserving 
the total zero momentum. Above the critical a\2 value the soliton with constituents of the 
smaller mass splits into two equal parts and the heavier component becomes a breather 
keeping its original place (third pair of panels). At the critical value of ayi both types of 
instability can be observed (third and fourth panel) simultaneously. 

Similar instability effects of a bright-bright and other soliton pairs have also been observed 
in Ref. 39]. The instability is explained by Kevrekidis at al. 40] as the appearance of a 
negative eigenvalue pair obtained by a linear stability analysis. It has also been shown that 
higher order bright-bright soliton pairs would exhibit instability irrespectively of the system 
parameters. Interestingly, it is possible for one of the higher-order bright solitons, i.e., for 
the one which has the stronger self-attraction, to recover its stability by collapsing into one 
bright soliton and expel the other component from its original position. 

Finally, it is important to note here, that the phases of the divided solitons are equal, 
therefore these particle-like wave-packets remain coherent with each other. Similarly to an 

15 



optical beam splitter where light is divided into two coherent beams, one could divide these 
matter waves and use one of them as a probe and the other one as a control packet. The 
probe packet could undergo transformations, while the control packet is left to evolve freely, 
thereby a phase-difference could build up between the two packets. If the two packets are 
brought together again, the phase-difference could cause interference pattern which allows 
one to quantify coherence. This may potentially be helpful for calibrational purposes as 
well. 



C. Dark-dark soliton pairs 

At the end we are examining the third possible combination of soliton pairs; a dark soliton 
is excited in each condensate. Interestingly this pairing can be stable even if the interaction 
inside each condensate is attractive (an, 022 < 0). The dark-dark solitons are described by 
the following formulas 



DD 



DD 



-iE® D tj exp 
exp ( —iE^ B t ) exp 



H 



H 2 



, &22C2 , ^ 

— O12C1 ) v z t 



3a) 
[23b) 



For our numerical investigations a Rb- Rb system has been chosen, with repulsive intra 



species interactions, an = 5.335 nm and 022 = 5.665 nm taken from [41]. This choice of 
coupling corresponds to scenario 4 of Table [1111 Tuning a\2 into the positive regime results in 
a stable dark-dark soliton pairs, see top two panels of Fig. HJ In the bright-bright case we have 
found that if the intra-species interactions are attractive and the inter-species interactions 
are repulsive, the solutions are unstable. We have tested numerically that this observation 
is valid in the dark-dark case as well, if the sign of the interactions is inverted. Furthermore, 
our numerical calculation provided an interesting phenomenon in the dynamics of these 
solitons. Preparing the dark solitons in the same way as before (i.e., by superimposing onto 
the ground state) but reversing the sign of ai2, resulted not only in losing the long-term 
stability, but also in the decay of a dark soliton in one of the components, and emerging as 
a secondary dark soliton in the other component (bottom two panels of Fig. HJ). 
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FIG. 4. (Color online) Evolution of a dark-dark soliton pair in a binary BEC composed of 87 Rb 
and 87 Rb atoms. Intraspecies scattering lengths are a\2 = 1.5 nm (612 = 3.62) (top panel), a\2 = 
— 1.5nm [b\2 = —3.62) (bottom panel). Other parameters are: an = 0.97 x 5.5 nm (fen = 13.04), 
022 = 1.03 x 5.5 nm (6 2 2 = 13.52), mi = 87 au, m 2 = 87 au, uj lt j_ = 2vr x 710 Hz, Ai = A 2 = 0.2, 
v = 0.2, Ni = 500, N 2 = 1600. The snapshots are taken at t = 50. 
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K 2 b u 


< b 2 2 


4d 






b\2 > b 2 2/K 


if 


K 2 bn 


< b 2 2 



TABLE III. Various domains of inter- and intra-atomic interaction strengths, bij, permitting the 
existence of a dark-dark soliton pair. 

IV. SUMMARY 

We have considered the existence and stability of soliton excitations in a two-component 
Bose-Einstein condensate both analytically and numerically. Our model allows the compo- 
nents to represent different elements (mi 7^ 777.2), but we also included those cases when two 
hyperfine states of the same element constitute the condensates. We excluded the possibil- 
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ity that these components can transmute into each other, i.e. the hyperfine states cannot 
be driven into each other. The dynamics of these condensates, within the mean-field zero- 
temperature approximation, are governed by the coupled Gross-Pitaevskii equations. We 
chose our presentation to be suitable for combining analytical results with earlier numerical 



investigations, e.g. 



26 



32 



42]. 



The occurrence of particle-like excitations together with conserved quantities are asso- 
ciated with the integrability of a nonlinear evolution field-equation, such as the coupled 
Gross-Pitaevskii equations. Using the results of a recent Painleve analysis of the coupled 
Gross-Pitaevskii equation {^fj, we showed how the system parameters determine the integra- 
bility of this system. However, for the CGP equations the studies so far restricted themselves 
for either equal coupling coefficients (fen = 612 = ^22) and/or equal masses. We note here 
that excitations with long lifetime may exist for a nonlinear evolution equation even when 
the integrability conditions are violated (43)], but these cases are possibly exceptional and 
do not represent the generic behaviour. 

We examined those ranges of system parameters which permit coupled soliton solutions 
of the governing equations fl5al -b). Here, we utilised the two well-known one-soliton solutions 
of the one-component nonlinear Schrodinger equation; the bright and dark solitons. These 
excitations then spatially modify the Thomas- Fermi ground state in our theoretical descrip- 
tion or by the appropriate ground state in our numerical simulation. We found analytically 
that each possible pair of these solitons have a range of system parameters where they are 
stable, and presented these ranges in tabulated form in Tabs. HEl and lHIl One can present 
these findings differently, if the intra-species coupling coefficients are held relatively fixed, 
and only the inter-species one is tuned. For example, if we assume bn and 622 to be positive 
(effective repulsive interaction) then the three tables [H HI1 and II III can be combined into one 
graph (see Fig. [5]). 

Moreover, we also examined how the stability of these soliton pairs is changing if the 
inter-species coupling coefficient is detuned from the value predicted by the Painlee-analysis 
via Eq. ([T|). It was shown, irrespectively of the pairing, that the stability is lost, although 
the pairs were not equally sensitive to the detuning, e.g. the motion of the bright-dark 
pair became erratic and preserved its periodicity only qualitatively after changing from 
5.5 nm to 5.4 nm. 

Two types of instability of static bright-bright soliton pairs composed of species with 
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FIG. 5. Ranges of the inter-species coupling coefficient, b\2, in which different types of soliton pairs 
can exist and are stable. The top graph is valid if K 2 bu < 622, while the lower graph is for the 
complementary case, i.e. K 2 bn > &22- The dotted line shows the result of the Painleve-test ([I]) in 
the case of harmonic trapping potentials. 



unequal masses have been observed. When the interspecies interaction lies below a critical 
value (0 < 612 < &12)' the static bright-bright pair evolves into a repulsive, momentum 
conserving, moving soliton pair. When the value of 612 > then one of the bright soliton 
(the constituent with smaller mass) splits into two equal portion of same phase while the 
other bright soliton becomes a breather. 

Well below the critical temperature of the Bose-Einstein condensate our description is 
expected to be adequate and the results could help experimentalist to modify the scatter- 
ing lengths via Feshbach resonance into a range where stable soliton pairs exist. As the 
temperature increases, however, the interaction with the thermal cloud becomes more and 
more important, and could not be neglected any more. We have not yet examined how 
the interplay between the condensates and the thermal clouds (for each component) would 
modify the dynamics. This needs further research beyond the mean field description. 
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